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Abstract 

Typically, the exact ground state energy of integrable models at finite volume can be computed 
using two main methods: the thermodynamic Bethe ansatz approach and the lattice discretization 
technique. For quantum sigma models (with non-ultra local Poisson structures) the bridge between 
these two approaches has only been done through numerical methods. We briefly review these two 
techniques on the example of the SU(2) principal chiral field model and derive a single integral 
equation based on the Faddeev-Reshetikhin discretization of the model. We show that this integral 
equation is equivalent to the single integral equation of Gromov, Kazakov and Vieira derived from 
the TBA approach. 



1 Introduction 



The spectrum of quantum integrable models in 1+1 dimensions can be computed exactly. In 
an infinite space (line) the spectrum is continuous and the relevant quantities are the dispersion 
relation e(p) of the fundamental particles and their two-body 5-matrix S(p,p'). 1 In a finite circle 
the spectrum of the system is quantized. When the circle length is very large, the momenta of the 
particles are constrained by a set of asymptotic Bethe equations 

N 

j" L l[S(p j ,p k ) = l (1) 

and the energy is given by the sum of the individual energies of the particles, E = Ylj=i e (Pj)- This 
treatment is not exact in the system size. If L is not very large, quantum corrections associated 
to virtual particles going around the space-time cylinder invalidate the asymptotic Bethe ansatz 
treatment. For relativistic systems, these corrections are typically or order e~ mL where m is the 
mass of the lightest particle [1, 2]. 

There are two main approaches to study the exact spectrum of quantum integrable theories at 
finite volume. One is based on a Wick rotation trick devised by Zamolodchikov [3]. The idea is to 
exchange time and space in the path integral. The finite circle length leads to finite temperature 
after Wick rotation. The infinite time over which we preform the path integral becomes an infinite 
space extent. In the Wick rotated picture the asymptotic Bethe ansatz description is therefore 
exact. This trick can be made rigorous for the computation of the ground state. Under some 
assumptions, it can be generalized to excited states [4, 5]. 

The second approach is based on integrable lattice discretizations of the quantum field theory. 
Frequently, we can interpret particles excitations of the integrable model as effective excitations 
over a non-trivial vacuum of a much simpler bare theory. Often, the bare theory is the continuum 
limit of a quantum spin chain. The effective particles of the field theory are the spinons of the 
quantum spin chain, i.e. the spin chain excitations around its anti-ferromagnetic vacuum. The 
bare theory is typically ultra local and admits an exact Bethe ansatz solution. The exact spectrum 
of the effective theory can be studied by carefully following the continuum limit of these exact 
equations. 

For many examples, it was understood how to match these two approaches analytically. However 
a systematic understanding of the connection of the two methods is still lacking. In particular, 
for quantum sigma models with non-ultra local Poisson structures the bridge between these two 
approaches has only been done through numerical methods. It is the purpose of this paper to 
preform a first analytical comparison of the two methods. For illustration we will consider the 
example of the exact ground state energy of the SU(2) principal chiral field [6, 7]. 

The exact vacuum energy 2 Eq of the SU(2) principal chiral field was recently considered by 
Gromov, Kazakov and Vieira (GKV) through the Wick rotation approach. The outcome of this 

1 We are oversimplifying since in general particles might have additional internal degrees of freedom and 
the S'-matrix is a matrix. The n-body S'-matrix factorizes for integrable models and therefore the main 
building block is the two-body 5-matrix. 

2 The authors considered not only the vacuum energy but also all the excited states however as explained 
above we will focus on the ground state energy in this work. 
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approach is a single integral equation given by 



g 2 = _ e imisinh(™) exp 2i j m 



Tfl f 

Eq = — — / ducosh(7ru) log 



(fl+) 2 - 1 (g-) 2 - 1 

\g+\2 _ 1 | 5 +|2 _ I 



(2) 



The star denotes convolution, the superscripts ± indicate shifts of the functions in the imaginary 
direction and Kq is the Kernel of the SU(2) S-matrix. These quantities are rigorously defined in 
the main sections of the text. 

The lattice discretization approach started with the seminal work of Faddeev and Reshetikhin 
[8]. In this work the authors proposed a discretization of the SU(2) principal chiral field, whose 
diagonalization is equivalent to the one of an inhomogeneous SU (2) spin s chain. This is quite an 
interesting proposal with several subtleties reviewed in section 3.1. The computation of the exact 
spectrum of the SU(2) principal chiral field can be done using this discretization. This was done 
in [9] 3 leading to a set of two coupled integral equations of the form 4 

loga = - mL cosh (mi) + -K * log(l + a ) - \ k q + * ^si 1 + «o) 

- - log(l + So) + s * log(l + o)(l + a) 

111 (3) 
log a = -Kq * log(l + a) - ^ K o + * lo &( 1 + a ) ~ 2 log ^ + a ) + s * lo s( 1 + a o)( 1 + a o) 



Tfl f 

E = —— / du cosh(iru) log (1 + a )(l + a ) 



Numerical analysis indicates that the Wick rotation approach and the lattice discretization approach 
agree [12, 13]. In this paper, we show that in fact this set of equations can be reduced to a single 
one. Moreover, we also show that this new single integral equation is equivalent to (2), once a 
proper change of variable is performed. 

We will also discuss what is the reason for the seemingly more complicated structure of (2) 
compared to the usual Destri-de Vega equation [14, 15]. 5 As explained in the last section this 
difference is directly related to the structure of the spin chain vacuum which is given by a condensate 
of bound-states known as Bethe strings. In the cases where DdV holds, for example for the SU(2) 
Gross-Neveau model, the vacuum distribution is typically given by a sea of real momenta particles. 

The structure of the paper is as follows: in section 2, we introduce the model and expose the 
necessary tools for the GKV method; the section 3 will be devoted to the origin of the lattice 
formulation of the model, and to the description of the procedure to obtain the corresponding 
integral equation; finally, in section 4, we show the equivalence of the two approaches and give 
some physical insight about the meaning of the obtained equations. 



3 Actually, the approach in [9] is based on the light-cone discretization proposed in [10], which turns out 

to be equivalent to the FR model [10, 11] 

4 The convolutions with Kq + appearing in these equations are understood in the principal part sense 
5 More precisely, in the usual DdV approach we typically obtain (2) without the denominator factors 

appearing in the logarithms of those equations. 
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2 Review of the Wick rotation approach 



2.1 The SU(2) Principal Chiral model 

The two dimensional SU(2) Principal Chiral model is defined by the following action 



S = "77^ / dtdxtr (g-%g) 2 , (4) 



1 

2~e 

with g £ SU{2). It has the global symmetry SU(2)l x SU(2)r which corresponds to the left and 
right multiplications by constant matrices go,g± € SU(2), i.e. g — > g§gg\- We consider the space 
direction to be a circle of length L. In infinite volume, the fundamental excitations are massive 
relativistic particles transforming in the bifundamental representation of SU(2) x SU{2). We use 
the rapidity 9 to parametrize the energy and momentum of these excitations as 

e(9)±p{9) = me ±7re . (5) 

The two body S'-matrix describing the scattering of these particles is completely fixed by Yang- 
Baxter, unitarity, crossing symmetry and absence of bound-states and is given by [16] 

S(e 1 ,e 2 ) = ^(ei + iF)®^(9I + i¥) where S (9) = ? || ~ fj ^ | +fj , (6) 

with 9 = 9 1 - 2 . 

Next we move to the large volume (Lm 3> 1) spectrum. We consider N particles in the circle. 
The simplest possible configuration for the SU(2)l x SU(2)r internal degrees of freedom is to have 
all the left spins and all the right spins pointing up. A generic configuration is obtained from this 
one by flipping Nl left spins and right spins. The large volume (Lm 3> 1) spectrum is obtained 
from the periodicity of the multi-particle wave function. This condition takes the form of a set of 
(nested) Bethe ansatz equations [17] that entangle the several degrees of freedom of the particles: 
the physical rapidities (O's) and the internal rapidities (u's and v's) (see figure 1) 

N N R ft a, M i Nl ft „ 4- ' 

e mL tt S 2 {e _ e ) rr g? ^±g TT °1 Vk + 2 = 1 (7) 

N Q i N R 

TT u 3 ~°k-2 = TT u i -U k -i (8) 
TT v 3 -Ok - 2 = XT Vj -Vk-l (Q) 

The asymptotic spectrum is simply the sum of the free dispersion relation of the particles 

N 

E = ^mcosh(7r/9 fc ) . (10) 

k=l 

For smaller volumes the asymptotic Bethe ansatz description breaks down and new techniques 
need to be developed. The simplest quantity we can consider in finite volume is the ground-state 
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Figure 1: Schematically, the SU(2) PCF describes massive particles 
with momenta parametrized by 9, and internal degrees of freedom 
corresponding to bifundamental representations of SU{2) x SU(2). 
Exciting these isotopic sectors is equivalent in this picture to flipping 
the corresponding spins. 



energy of the theory, or Casimir energy. The most common procedure to compute this quantity is 
by means of the Thermodynamic Bethe Ansatz (TBA), which was originally developed to compute 
finite temperature thermodynamics of integrable models, i.e. its free energy. If we perform a double 
Wick rotation and simultaneously interchange the role of time and space variables we realize that 
the free energy per unit length at temperature 1/L actually coincides with the ground-state energy 
at finite volume L [3]. 6 

Once we consider the system at finite temperature and infinite length we excite an infinite 
number of all possible excitations, N, Nl,Nr — > oo. In this limit, if some rapidity Uj has an 
imaginary part, the left hand side of (8) diverges (or decays to zero). The right hand side must 
have the same behaviour. Hence, there ought to be another rapidity separated by i relatively to 
Uj. Therefore, it is admissible that all Bethe roots organize in n-strings, i.e., a set of n Bethe roots 
with the imaginary part separated by i (this is known as the string hypothesis [18]). Obviously, the 
same argument could be applied to the other wing associated to v rapidities. At finite temperature 
we should expect infinitely many strings of all possible sizes. The proper variables to describe these 
degrees of freedom are therefore densities of Bethe roots and holes [19] associated to each of the 
strings and physical rapidities. The asymptotic Bethe equations ought to be written using these 
densities. We introduce the densities p n and p n of Bethe roots and holes, reserving the index n = 
for densities of 9's, n > for the densities of n-strings of it's, and n < for strings of v's. 



6 The double Wick rotation does not change the theory due to Lorentz invariancc. For non-Lorentz 
invariant theories the Wick rotation trick can often still be applied but it would relate the free energy of one 
theory with the ground-state energy of another theory. 



The Bethe equations (7), (8), (9) become [12] 




(11) 
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where * stands for the usual convolution 

+00 

/*<?= J d0'f(0-0')g(0') , (12) 

—00 

The kernels K nm are the derivative of the logarithm of the S-matrices between the several excitations 
{u strings, v strings and physical rapidities), see [12] for details. For instance, 

-K 0t0 {6) = K (0) = ^^log^(e) . (13) 

To find the ground-state at finite temperature (1/L) and infinite volume, we must minimize the 
free energy. The outcome is a set of equations constraining the densities, known as TBA equations 

logy n + 5 n fimL cosh ttu = s * log(l + Y n+1 )(l + Y n -{) , (14) 

where s = 2 cosh mi • These equations are written in terms of some thermodynamic functions Y n . 
In the ground state, they are simply related to the densities by logy n = p n /Pn for n 7^ and 
log Yo = po/po- Moreover, the exact finite volume energy of the system is given by 

777 f 

E (L) = - — / ducosh(7ra)log(l + Y ) . (15) 

Equivalently, if we assume that \ogY n (u) has no singularities in the physical strip \Im(u)\ < i, for 
any n, then the TBA equations are equivalent to the so called Y-system 

y+y- = (1 + y n+ i)(i + y n _o , (16) 

with the asymptotic "boundary conditions" Y n (u) oc e~ mLS "'° cosh7rn for large u. Throughout this 
paper, we will often use the notation 

/± = f(u ± l -) , /±± = f(u ± i) , = f{u ± |) . (17) 

2.2 The Wick rotation method 

Using the TBA approach we obtained the value of the ground state energy of the system at finite 
volume. However, the solution involves an infinite number of coupled equations. The GKV method 
[12] allows one to collapse these equations into a single integral equation. 7 This method in strongly 
based on the Hirota equation 

T+T- - T„_ 1 T n+1 = $N$hn] ) (18) 
which is related to the y-system (16) by the change of variables 

_ r re+ ir ra -i 



furthermore, it can be used to study excited states of the theory in a systematic fashion. 
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The Y-system is actually a gauge invariant version of the Hirota equation, as can be seen by noting 
that it remains unchanged by the transformations 



T k -> g [k] g [ - k] T k 
<3? ->■ g + g~<& 
$ ->■ g + #~<I>, 



(20) 



where g and g are arbitrary functions. To ensure reality, g is the conjugate of g and is the 
conjugate of <I>. 

The crucial observation is that the general solution to Hirota equation (18) is known! It involves 
a two by two determinant with four unknown functions, see e.g. [20]. Requiring a real solution we 
reduce this number to two functions (and their complex conjugates). This means that effectively 
we only need to find these two functions to solve (18) completely. Furthermore, as described above, 
there is a huge gauge symmetry of Hirota which relates different solutions to Hirota equation 
yielding the same Y- functions. Using this gauge symmetry we can render one of the two remaining 
functions trivial 8 . Thus, at the end of the day there is one function left unfixed. We can trade 
this function for another interesting unknown as we now describe. For the vacuum, the Y-functions 
for the two wings must be symmetric, Y k = Y_ k . This does not necessarily imply that T k = 
but it does implies that T k and T_ k are related by a gauge transformation g(u). The knowledge of 
this function is equivalent to the knowledge of the last unknown function in the general solution of 
Hirota. At the end of the day, the infinite set of integral equations for the Y n functions become a 
single integral equation for the gauge function g. It reads [12] 



imL sinh(7nr) 



exp 2i Im 



K Q *log 



(g + Y 



i 



7+12 



(21) 



The ground-state energy can be expressed in terms of this gauge function g, as 



E (L) 



m 
~2 



du cosh (7™) log 



(g+f - 1 {g- Y 



Iff 



1 b+| 2 -l 



(22) 



Therefore, the problem of computation of the exact ground-state energy at finite volume is now 
solved by a single equation. 

This method can in principle be applied to many integrable theories. For instance, the finite 
volume ground-state of the model SU(2) Chiral Gross-Neveu [21] was solved by this method, 
resulting in the well-known Destri-de Vega (DdV) equation [14, 15] . This equation turns out to be 
exactly the same as (21) if we simply drop the denominator |<? + | 2 — 1 inside the logarithm to be 
left with 
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imL sinh (ttx) 



exp (ilm K * log ((5^ 



!))• 



(23) 



The ground-state energy can again be expressed in terms of g. Again, the result turns out to be 
exactly the same as in (22) once we drop the denominator \g + \ 2 — 1 in this expression. In the last 
section, we will give a physical explanation for this distinction. 



^constant or a simple polynomial, depending on the situation, see [12] for details. 
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3 Lattice approach 



3.1 Faddeev-Reshetikhin model 

An alternative method for the quantization of the SU(2) PCF was proposed in [8] by Faddeev 
and Reshetikhin (FR). It is based on an integrable discretization of the model. In this approach 
the SU(2) principal chiral field particles are the effective excitations around the anti-ferromagnetic 
vacuum of a quantum spin chain. 

The direct quantization by the quantum inverse scattering method [22, 23, 24] of the action (4) 
turns out to be problematic due to the structure of Poisson brackets. Let us recall why and review 
what was the approach of FR to circumvent this point. We introduce the currents = d^gg^ 1 = 

L^t" 1 , where t a are the generators of SU(2) 9 . By definition, they satisfy the flatness condition 

d,L a v -dvLl + ie abc LlLl = U. (24) 

We can rewrite the action (4) in terms of these currents, introducing a Lagrange multiplier field 
2e^-K a (x) to fix the constraint (24) 

S = "2i§ E / dtdx (( L S) 2 + (LI) 2 + 2eg 7 r a (d 1 L§ - d Lf + i^Lg)) . (25) 

The current Lq is not dynamical in the sense that it appears with no time derivatives. The equation 
of motion for this current reads 

L = e ze ir L 1 — e^aiir . (Zb) 

We can now compute the Poisson brackets of the current elements, such as {Lq(x), L\(y)}. We 
simply need to notice that 7r a is the conjugated momentum to Lf so that the equal time Poisson 
bracket between these two objects is simply {Lf(x), ir b (y)} = S ab S(x — y). Therefore 

{L », L\{y)} = elie ahc L\{x)5{x - y) - e 2 5 ab 5'(x - y) . (27) 

The derivative of the delta function term rules out the application of the quantum inverse scattering 
method [22, 23, 24], see [25, 26, 27, 28] for further discussions of these nondocal terms. 

To overcome the difficulties generated by this term, FR proposed two inspired steps [8] 

• First, the current Poisson brackets are modified. The derivative of delta function terms are 
dropped. 

• Next, a new Hamiltonian is proposed. The new Hamiltonian yields the original equations of 
motion using the new Poisson brackets. 

Ultimately, the validity of these two seemingly dangerous steps is justified by the final results which 
match perfectly with all other approaches. The great advantage of the new Hamiltonian is that 
is allows for an integrable discretization. Physically, this Hamiltonian describes two interacting 
magnets. Once we discretize it, the diagonalization of the lattice model is equivalent to the one of 

9 These generators satisfy [t a ,t b ] = ie abc t c and tr t a t b = 5 ab 



an inhomogeneous SU{2) spin-s chain with N sites (N is even) and inhomogeneity parameter a. 
These parameters need to be sent to infinity in a precise way which we review below. 

For any finite s, N and a the bare spin chain model is described by an exact set of Bethe 
equations. The idea of the lattice approach to the finite size spectrum of the SU (2) principal chiral 
field is to use these exact equations and follow all the continuum limits carefully to derive all finite 
size corrections to the effective model. 

The Bethe equations for the inhomogeneous SU(2) spin-s chain with N sites and M Bethe roots 
[8] read 

— / km 
/ Uj + a + is\ 2 / Uj - a + is\ 2 _ -i-r Uj - u k + i 

\Uj + a — is J \uj — a — is J uj — u& — i' 

where a plays the role of the inhomogeneity. These are the bare Bethe equations describing the 
magnonic excitations over the ferromagnetic vacuum. To make the bridge with the SU(2) PCF we 
need to consider the solutions corresponding to fluctuations around the anti- ferromagnetic vacuum. 

In the large N limit, the anti- ferromagnetic vacuum of this chain is composed by a Dirac sea 
of N/2 2s-strings of Bethe roots (M = Ns) [29, 30]. The claim is that the excitations around this 
vacuum are the SU(2) PCF excitations. The limit is taken as follows 



• The physical length of the model 10 , mL, is given by 

mL = 2Ne-™ . (29) 

We should take N, a — > oo keeping this quantity fixed. In this limit, the holes in this vacuum 
become massive excitations with relativistic dispersion relation. 

• The scattering matrix for these holes becomes the scattering matrix of the SU{2) PCF when 
we take furthermore the large spin limit, 11 

s -> oo . (30) 



As mentioned above, the physical rapidities of the SU{2) PCF are identified with the hole excitations 
of the anti- ferromagnetic vacuum that we denote by u 2s (see appendix A) , 

9 o u 2s . (31) 

To describe the physical excitations of the SU(2) principal chiral field we also need the rapidities u 
and v introduced in (7)-(9). These turn out to correspond to string excitations of sizes 2s + 1 and 
2s — 1 respectively, whose centers we denote by u 2s+1 and m 2s_1 , 

u 4* u 2s+1 (32) 
v o u 23 " 1 . (33) 

More precisely, the equations for u 2s , u 2s_1 and u 2s+1 can be derived from (28). For large length Lm 
they coincide with (7)-(9) provided we identity u 2s , n 2s_1 and u 2s+1 with 9, v and u, respectively 

10 This is the length of the space circle measured in unit of inverse infinite volume mass gap 
11 It would be very interesting to explore what the finite s model would describe and what the connection 
between the lattice and the Wick rotation approaches in that case. 
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(see also [31, 32]) 12 . This is reviewed in appendix A using the DdV approach. It is the first major 
evidence for the correctness of the FR discretization! We observe however that the FR model 
generates only a part of the Hilbert space of the SU(2) PCF. As showed in [8], the number of holes 
is always even and only singlet states in one isotopic sector can be excited. 

strings holes strings 

2s-3 2s-2 2s-l 2s 2s+l 2s+2 2s+3 ••• s — ► oo 

t 

T 

v e u 

-3 -2 -1 1 2 3 

Figure 2: In infinite volume, the degrees of freedom of the spin 
chain have a correspondence with those of the SU(2) PCF. However, 
a counting of the states of the spin chain shows that only singlets in 
one of the SU (2) sectors are generated 



3.2 Finite Volume NLIE 

There are several techniques to study the finite size/temperature problem in quantum spin chains: 
the Destri-de Vega approach [14, 15], the nonlinear integral equation (NLIE) method suggested in 
[33, 34], the fusion procedure based [35, 36, 37, 38, 39], and an hybrid formulation based on both 
fusion procedure and NLIE [40]. When the spins of the chain are in higher representations there 
is an important new feature of the anti-ferromagnetic vacuum: it is made out of a sea of bound- 
states or Bethe strings. 13 This is the case for example for the FR model presented above. This 
characteristic of the models renders the hybrid formulation more suitable [40]; the DdV approach 
is not yet developed to comprise this case. 

The bare Bethe equations (28) follow from the diagonalization of a transfer matrix T2s [8]- The 
exact (bare) energy of the system is also given in terms of the transfer matrix. Interestingly, this 
transfer matrix is only one in a family of fused transfer matrices T p with p = 1,2,.... All these 
transfer matrices are related by simple functional relations. Together with the analytic properties 
of the transfer matrices, these functional relations have a unique solution. More precisely, the fusion 
procedure based method [35, 36, 37, 38, 39] converts these equations into integral equations of TBA 
type. Typically we end up with an infinite set of integral equations. Often these can be simplified 
into a finite set of equations [40]. 

In this section, we will start by reviewing how to obtain an infinite set of coupled equations of 
TBA type, using the fusion procedure. We will then derive a single non linear integral equation 

12 Rigorously, this correspondence is valid when u and v are real. 

13 In fact, the Bethe roots condense into strings only in the large volume limit. Otherwise deviations to 
this behaviour appear. 
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that allows to compute the finite volume ground-state energy of the SU(2) PCF. Our method is 
based on the proposals of [40, 9]. In [9], it was obtained a set of two coupled integral equations. 
But, as we will see, this can be reduced to a single one equation. 

The starting point of the fusion procedure is the introduction of a family of fused transfer 
matrices defined by 

■ Cj) 



T^E A m ( 34 ) 



i=i 



where 

p- 1 l ~ 2 nfa+iW-p-i] 



m=l— 1 n=0 



In this expression, 4>{u) = (u — a) "2" (u + a) "2" and Q(u) = Y\j(u — Uj) is the Baxter polynomial. 
The Bethe equations (28) ensure that functions Tj have no poles. Next we introduce another set 
of functions, the Y-functions. They are given by 



T i-i T i+i 



Y,- = —. , J ~ L , ~ • (36) 

Y[l=i ( t )[2p ~ j+2s]( f )[ ~ 2p+i ~ 2s] 

From (34) the following alternative definition also follows: 

T+TT 

1 + (u) = — : 3 -^- . (37) 

By definition, they obey the Y-system functional relations 

Yi-Yj = (l + Y i _i)(l + Y i+ i). (38) 

Recall that exactly the same equations arise in the TBA context, see (16). Of course, the derivation 
of the Y-system equations in both approaches is radically different. 

Once supplemented with simple analytic properties of the Y-functions, the Y-system equations 
can be converted into integral equations for Y n which admit a unique solution! The analytic 
properties of the Y-functions follow from those of Tj, which in turn depend on Baxter polynomials 
Q. In the large N limit, the zeroes of Q are given in terms of 2s-strings, and, for finite N, we expect 
the deviation of the imaginary part of these zeroes to be small (less than | according to [41]). 

The determination of the zeroes of Tj in the ground state needs some care. The zeroes can 
appear in two different ways. First, there are zeroes coming from the factorization of some 0's 
that are common in all terms of the expression for T,-, (34). For clarity, we call these zeroes the 
trivial zeroes. On the other hand, there are zeroes that appear due to the Baxter polynomials in 
the expression of Tj . 

The trivial zeroes appear only when j > 2s. When j < 2s there are no common ^'s in all terms 
of (34). So, if we let j = 2s + k (k > 0), then the trivial zeroes are the zeros of the common factor 



fe-i 

j-j" ^[fc-l-2m] 



(39) 

m=0 



The remaining zeroes can be determined by considering the large iV behaviour of T,- . Such kind 
of analysis is nicely discussed in [32] which we will follow closely. For finite N we expect a small 
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deviation from this behaviour as stated before. In the large N limit, either the first or the last 
term of (34) is dominant, depending on whether the imaginary part of the argument is negative or 
positive, respectively. In this way, we can obtain the approximate positions of the complex zeroes. 
Explicitly, if Im u > 0, 

TiW~|^n^ 2n+2s+i_i] (40) 



n=0 



and if Im u < 



^ n=0 

We now consider the ground-state configuration at large N, where the Bethe roots organize in 
2s-strings. Therefore, for j < 2s, some of the zeroes of the Baxter polynomial in the numerator 
cancel with the same zeroes of the Baxter polynomial in the denominator. Therefore, the position 
of the imaginary part of the zeros of Tj in the large N limit is given by 

s + i-l) with I = 0, - 1, forj<2s (42) 

a + i~l) with / = 0, ...,2s - 1, for j > 2s . (43) 

With these considerations we found Nj zeroes of Tj, for any j. Indeed, for j < 2s, all zeroes come 
from the Baxter polynomials. Taking into account (42) and the fact that in ground-state we have 
Y 2s-strings, the number of complex zeroes of Tj<2s is y x 2j = Nj. For j > 2s, we have N(j — 2s) 
trivial zeroes according to (39), and 2sN complex zeroes due to (42). This gives again a total of 
Nj zeroes. 

Moreover, a comparison of the number of zeroes in both sides of (34) shows that in the ground- 
state we do not expect extra zeroes on the real axis. Indeed, let rih be the number of real zeroes 
of Tj excluding trivial zeroes, so that its total number is rih + Nj, for any j. From (34), we see 
that the right hand side is the ratio of two polynomials of degree 2{j + l)Ns + Nj and 2(j + l)Ns. 
Furthermore, due to Bethe equations (28) we know that the zeroes in the denominator are cancelled 
by the zeroes in the numerator. Therefore, the right hand side of (34) is a polynomial of degree Nj 
due to the factors of 4> in each term A] , . Equating the number of zeroes of both sides gives that 
n/j = and we conclude that in ground-state there are no zeroes on the real axis. 

Using this information we can determine the region where Tj are analytic, non-zero with con- 
stant asymptotic behaviour (ANZC). The properties of Yj can be extracted from those of Tj using 
the expression (36). The functions Yj do not have zeroes nor poles in the strip |Im u\ < §, except 

N 
2 



for Y2s, which has ^--fold zeroes at the positions ±a. This follows from the observation that all 



trivial zeroes of Tj in the strip |Im u\ < | cancel with similar zeroes of the denominator in (36), 
except for the case Y2 S . Recall that a is going to be sent to infinity in a particular scaling limit but 
it is crucial to keep these zeroes at this point; as anticipated above they will be responsible for the 
dynamical mass generation of the effective excitations. All these ANZC regions are summarized in 
the figure 3. 

We will derive integral equations for Yj from the Y-system (38). Let us define a new function 
Y2s = which is now ANZC in the strip |Im u\ < |. This function satisfies the equation 

V^ ™^> . ,44) 
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Tj<2 s (!i) and Y m ,{tl) T ]>2s (u) Y 2 .(u) 
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Figure 3: The regions ANZC of the functions Tj and Yj are repre- 
sented in grey, the black dots represent ^-fold isolated zeroes and a 
is the inhomogeneity parameter. 



We now consider the shift operator s = 2 C osh ttu w ^ich has the following property: for an analytic 
and bounded function f(u) in the strip |Im u\ < ^, we have 

s*[f + + r] = f. (45) 

We take the logarithm and apply the shift operator in both sides of (38) for j ^ 2s to get 

f logYi= s*log(l + Y 2 ) 

{ 46 
[logY,= a*bg(H-Y i _i)(l + Y J - + i), j^2s. 

Similarly, doing the same to equation (44) we get 

logY 2s = log^- S *log0 + <T + S *(l + Y j _ 1 )(l+Y J+1 ). (47) 

When we consider the continuum limit in which a — > oo and N — s> oo, with the quantity mL = 
2Ne~ na held fixed, the terms involving <^'s become —mL cosh (ttu). Finally the equation is written 

as 

logY 2s = -mL cosh(^u) + s * (1 + Y i _i)(l + Y j+1 ) . (48) 

The ground-state energy can be expressed in terms of Y 2s , by the following considerations. 
According to [8], the energy of a state is given by 

E = 1 log Q (~ a + ^ Q{a ~ ^ , (49) 
id Q(—a — is)Q(a + is) , 

where we have introduced the lattice spacing 5 to anticipate the continuum limit. Using the 
expression (34), it is possible to express the Baxter polynomials in terms of T 2s and <f>. The reason 
is that when T 2s is evaluated at ±o ± i/2 almost all terms vanish due to </>(±a) = 0. It remains a 
single one which involves precisely the Baxter polynomials appearing in (49). Then, we get 

E= l^^SL + ^ ^!±^ . (50, 
6 T+(-a) 5 ^ 4>{a + i{n + I)) 

The term involving fts is present in all states and gives an infinite constant contribution when the 
continuum limit is taken. Therefore, we will ignore it. Focusing only on the first term of (50) we 
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can rewrite it using the previous analytic properties of T 2s in the ground-state. In particular, with 
the equation (36) we can write 

= s * log(l + Y 2s ) + terms involving 0's . (51) 

Plugging this result in (50) and ignoring once more the divergent terms involving 0's, we finally 
obtain the ground-state energy in the continuum limit 

+oo 

£ = -y / dncosh(vrn)log(l + Y 2s {u)) . (52) 

— oo 

As before, we took the continuum limit by letting 5 — > and a — > oo with the quantity niL = 
2Ne~ 7va held fixed. Thus, we derived a system of equations resembling the TBA. Indeed, when 
s — >• oo the identification Y 2s+J - = Yj leads to a match of both approaches. 

We want now to truncate this system of equations using only functions which are natural from 
the lattice discretization point of view. The first step in this direction was given by Suzuki [40], by 
the introduction of the following function 

a (u) = . (53) 

A 2s+l\ U + 2> 

This function (and its complex conjugated) allows to truncate the number of equations in the 
Y-system as we will see. It can be written in the following way 

a (u) = T 2s . 1 q[ _ 2s] n 2,-i 0[an+2] - ( 54 ) 



A consequence of the definition is 

n':=oV [2m+2 ] q [ ~ 2s] 



T + q[2s\ 

1 + °oW= 2-1 2 »Z_^ (55) 



and also of great importance is the relation 

(l + ao)(l + oo) = 1 + Y 2s , (56) 

which is again a simple consequence of the definition of ao and of Y 2s . The expression (56) truncates 
the number of equations in the Y-system. Thus, one can think of the functions ao (and its complex 
conjugate) as a re-summation of the contribution of the Y functions with index larger than 2s — 1. 
Moreover, we see that the ground-state energy can also be expressed in terms of ao as 



En 

u 2 



m 



j du cosh(7ru) log (1 + ao)(l + ao) . (57) 



From the definition, one can determine the ANZC regions of the functions ao and 1 + ao which 
will be useful for what follows. They are summarized in the figure 4. 

With these analytic properties, we are in condition to perform the trick of [40]. Consider the 

3 

r 



function T 2s which is ANZC in the region |Im u\ < §, then by Cauchy's theorem 
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Figure 4: The regions ANZC of the functions 1 + ao(u) and ao(u) 
are represented in grey, the black dots represent ^-fold isolated zeroes 
and a is the inhomogeneity parameter. 



Figure 5: The contour T encircles an ANZC region of the function 
T 2s (z) 



<{dz dlOg J M e-^ = 0, (58) 
Jr dz 

where T encircles the real axis counter-clockwise, in the upper plane with z = | and in the lower 
plan with z = — \ (see figure 5). 

The relation (58) allows one to express the Baxter polynomials Q in terms of 1 + ao and 1 + clq 
using the expression (55). As explained in appendix B, we obtain 

logYi = s*log(l + Y 2 ) 

logYj- = S *log(l+Y i _ 1 )(l+Y ?+1 ), i = l,...,2s-2 
< logY 2s _i = s*log(l + Y 2s _ 2 ) + s*log(l + ao)(l + ao) (59) 
loga = —mL cosh 7ru + \Kq * log(l + ao) — \K^ + * log(l + ao) — \ log(l + ao) 
+s*log(l + Y 2s _i), 

where we have taken as before the continuum limit when computing the terms involving the several 
(/>'s. Also, the convolutions with Kq + are understood in the principal part sense. Note that at this 
point we did not take any limit on the spin s; it can be arbitrary so far. 

As observed by Hegedus [9], the limit s — > oo leads to an extra symmetry, due to the fact that 
the Y-system becomes symmetric with respect to the massive node, Y 2s , 

Y 2s+1 = Y 2s _x . (60) 
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This observation can be used to reduce the previous set of equations to a single equation for a . 
To derive this we observe that 



s * log(l + Y 2s -i) 



1 

2 
1 
2 

- (log((l + a )(l + a ) - 1) + mLcosh(7ra)) 



s*log(l + Y 2s+1 )(l + Y 2s _ 1 ; 
(log Y 2s + mL cosh(-/rti)) 



(61) 
(62) 
(63) 



The first equality uses Y 2s _i = Y 2s+ i which is the symmetry restored as s — > oo as described 
above. In the second equality, we used (48). The last equality follows from (56). Therefore the 
system of equations (59) can be simply rewritten as 



a = e -^ cosh ^exp 



X -K Q * log(l + oq) - ^K++ * log(l + d ) + - log((l + a ) - (1 + ao)" 1 



(64) 



4 Equivalence of both approaches 

As reviewed in the previous sections, the Wick rotation approach and the lattice discretization 
approach to the finite size problem of the SU(2) PCF lead to two seemingly distinct set of integral 
equations, (2) and (3). In [12, 13] numerical evidence was given in favour of their equivalence. The 
purpose of this section is to provide an analytic prove of their interchangeability. 

Comparing the formulas for the energy of the ground state given in both approaches, (2) and 
(57), we are inspired to try 

(g+) 2 - 1 

1 + a = TTT? T ' ( 65 ) 

so that the expressions for the ground state energy take the same functional form. Indeed, with 
this simple change of variables we can prove the equivalence of the two approaches. First let us 
re-write the GKV equation (2) by shifting the contour of integration as 

G? + ) 2 = -e-— exp (k * log ^1 - K++ * log - log |^|) . (66) 

This is actually the form of the equation which is most suited for numerics [12]. Convolutions with 
Kq + are understood as principal value integration. Under the identification (65), the equation (64) 
coincides precisely with (66)! 

We observe that the approach followed in [9], resulting in equations (3), is equivalent to ours. 
Indeed, the equations (3) can be reduced to a single one. According to [9], the functions a and a 
satisfy 

(l + a)(l + S) = l + Y 2a +i. (67) 
Using the trick described in (61-63) we obtain precisely (64). 
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5 Conclusion and discussion 



We have shown analytically the equivalence between the approaches studied in this paper: the Wick 
rotation and the integrable discretization based methods. Indeed, the change of variable (65) leads 
to a perfect match between equations (2) and (64). This also confirms the numerical agreement 
presented in [12, 13]. 

The GKV equation for the SU(2) PCF (21) has a striking similarity with the typical DdV 
equation (23), with an important difference: the denominator |<7 + | 2 — 1 appearing in the convolutions 
of (2) is absent in (23). The DdV approach converts the bare Bethe equations of a lattice model 
into an exact integral equation for a counting function (see appendix A). This integral equation 
can be used to compute the exact anti-ferromagnetic ground-state energy for any number of sites, 
N. When the continuum limit of the spin chain is considered, this approach gives the finite volume 
corrections of the corresponding field theory. Typically g 2 appearing in (2) and (23) can be identified 
with the counting function in the continuum limit. Unfortunately, this method is designed to study 
only lattice models with anti-ferromagnetic vacuum made of real Bethe roots. 

When the vacuum is a sea of complex Bethe roots, the DdV technique is not yet developed to 
comprise this case. Generally, the study of spin chains at very large A is made using the string 
hypothesis. This essentially means that all the complex Bethe roots organize themselves into sets 
with the same real part (the center of the string) and imaginary parts separated by i. When the 
number of sites is not very large, it is assumed that the corrections are exponentially suppressed in 
A. 

In principle, we could fuse the bare FR Bethe equations (28) to obtain some equations for 
the centers of the strings. Since the centers are real, the DdV is suitable. More importantly, the 
corrections given by DdV are exponentially suppressed in mL ~ Ae^ 7 ™ <C A, so we could expect 
that the string hypothesis would be justified. Such approach is done in appendix A, but the result 
is similar to (23), i.e., the GKV equation (2) without the denominator \g + \ 2 — 1. 

The explanation for this mismatch is that the deviation to the string hypothesis is not negligible, 
because the number of Bethe roots M = As also goes to infinite. This is a somehow subtle point: 
when we have A — > oo and a finite number of roots which form strings, these strings are extremely 
vertical. More precisely, the deviations from the string hypothesis are exponentially suppressed in 
A S> Ae _7ra . However, when the number of Bethe roots also scales with A" the suppression becomes 
much smaller and the string hypothesis cannot be used, see also [41]. In other words, we see that 
the denominators in the GKV equation for the SU(2) PCF have a very nice physical meaning: 
they are a manifestation of the string hypothesis deviation in the vacuum of the underlying lattice 
model 14 . 

Unfortunately, the FR model does not generate the whole Hilbert space of the SU(2) PCF 
[11, 8]: the number of holes is always even and only the singlets in one isotopic sector can be 
excited. Nevertheless, a modification of the NLIE (3) to include the one particle excited states was 
proposed and numerically tested with success [42, 13]. It would be interesting to check if it still 
agrees with the GKV approach which is claimed to be valid for all excited states. We also hope that 

14 In the cases where the DdV method is suitable, as in the SU{2) Chiral Gross-Neveu, the exact DdV 
equation is recovered by the Wick rotation approach. It is well known that the SU (2) Chiral Gross-Neveu 
model has a lattice formulation through the simple XXX i inhomogeneous spin chain (see for instance [32]). 
The anti-ferromagnetic vacuum of this spin chain is made of real Bethe roots and, consequently, the DdV 
can be directly applied to the bare Bethe equations. 



17 



the comparison between the two approaches made in this paper can help to determine a perfect 
discretization of the model, comprising all the excited states. 

One of the main motivations for the study of the present paper came from the study of the finite 
size spectrum of planar AdS / CFT describing strings in AdS$ x S 5 and N = 4 supersymmetric Yang- 
Mills theory. The asymptotic spectrum is given by Beisert-Eden-Staudacher equations [43, 44], the 
finite size Y-system was proposed in [45] and the corresponding TBA equations were presented in 
[46, 47, 48]. These three results are the analogue of equations (7-9), (16) and (14) for the SU(2) 
PCF. The analogue of (2), i.e., a finite number of equations describing the full planar spectrum of 
AdS/CFT, is still not known despite many interesting recent advances [49, 50, 51]. It would be 
interesting to also try to look for these finite set of equations from the lattice discretization point of 
view. Furthermore, if such DdV like equations are found from the Wick rotation method it would 
be interesting to understand their physical meaning. Will they contain denominator terms like (2)? 
If so, do they signal some intricate structure of a lattice vacuum as in the case of the SU(2) PCF? 
As a warm up, and to gain some more physical intuition about the form of the final set of non-linear 
equations, we could generalize the lattice method to the SU(N) PCF which was considered using 
the Wick rotation method in [52, 32, 31]. 
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A DdV applied to the fused FR Bethe equations 

In this appendix, we will analyse the FR Bethe ansatz equations (28) both in infinite and finite size, 
using the Destri-de Vega approach. First, we want to study the large volume solutions of the FR 
Bethe ansatz equations (28) corresponding to fluctuations around the anti-ferromagnetic vacuum. 
We expect these effective excitations to have the same dynamics as the excitations of the SU (2) 
PCF. To accomplish this, we will use the DdV approach. This technique is often applied to the 
computation of finite volume spectrum of quantum field theories with a discretization admitting 
an anti-ferromagnetic vacuum made of real Bethe roots. As we saw, our spin chain is a sea of 
complex Bethe roots which condense into 2s-strings, when the number of sites N is very large. In 
this limit, we can then fuse the bare Bethe equations, which amounts to multiplying the equations 
associated to the same string. In the end, we will get equations for the centers of the strings. We 
will apply the DdV technique to these fused equations, and ignore the finite volume corrections. In 
the continuum limit, we expect therefore to obtain the large volume description of the SU(2) PCF, 
i.e., the Asymptotic Bethe ansatz equations. 

The finite size study of these equations by the DdV approach has an instructive premeditation. 
In the string hypothesis, it is assumed that the corrections to it are exponentially suppressed in N. 
As we will see below, the DdV approach applied to the fused FR bare Bethe equations, gives, in 
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the continuum limit, corrections exponentially suppressed in ~ Ne~ na <C N. Therefore, a priori, 
we could expect that the string hypothesis deviation should be subleading compared to the DdV 
finite volume corrections. However, it will turn out that the string hypothesis deviation is not 
negligible, because the number of Bethe roots M = Ns also goes to infinity. In other words, the 
string hypothesis is no longer exponentially precise in N when the number of strings also grows 
with N. The effect of the curvature of the strings will be then manifest, when we compare the 
result with the correct GKV equation (21). 

In the large N limit, the fused FR Bethe ansatz equations (28) are obtained by multiplying the 
Bethe roots belonging to the same string. After taking the logarithm, we obtain 

AT 00 £™ 

y {e n ,2.(u? ~ a ) + d n,2s(u] + «)) = 2 ^Q] + E E S ">™K " u k) > ( 68 ) 

m=l k=l 

where u™ is a center of a string of size m, £ m denotes the number of strings of length m and is 
an integer or half-integer (depending on the particular configuration considered). 

The functions # n ,m an d E n ,m have the following expressions 




The integrals are understood in the principal part sense, due to the singularity at the origin. 

For simplicity, we will consider a state with M strings of size 2s, Nr strings of size 2s + 1 and 
Nl strings of size 2s — 1, with Q" being integer. Anticipating a result, we will denote uff +1 by Uf~, 
u k S ~ 1 by v k and the holes by Ok- We recall that holes correspond to the solutions of equation 
(68), distinct from the Bethe roots. We introduce the simpler notation 



V h {u) 


= 2s ,2s(u - a) + 2s ,2s(u + a) 


V u {u) 


= 2s +i,2s(u - a) + 2s+ i,2 S (u + a) 


V v (u) 


= 2s ~i )2s (u - a) + 02s-i,2s(« + a) 


H{u) 


= ^2s,2s(u) 


U(u) 


= ^2s,2s+l(u) = E 2s+ i :2s (u) 


U u (u) 


= S2s+l,2s+l( tl ) 


U v (u) 


= E2s+l,2s-lW 


V(u) 


= S 23) 2a-l(«) • 



Following the standard DdV approach, we define the counting function Z(u) as 

N M N R N L 

Z(u)= jV h (u)-J2nu-ut s )-J2^-^)-^2v(u-v k ). (71) 

k=l k=l k=l 



19 





k 1 ?! 

ie 








< 








-#-o — • — •— o — •-• 


► 




-ie 7 



Figure 6: The DdV contour 7 encircles an analytic region (with non 
divergent asymptotic behaviour) of the function B,2 S ,2s( v )- The black 
dots represent the centers of 2s-strings and the white dots represent 
the holes. 



The counting function has the property of Z(u k ) = 2irQ k ]s , where u k is a center of a 2s-string 
or a hole. Essentially, the DdV technique allows to convert a sum over the centers of 2s-strings 
into a sum over the holes in an exact way. In other words, using it, we can obtain the effective 
Bethe equations describing the excitations around the anti-ferromagnetic vacuum. Hence, it is an 
appropriate way of getting information about the dynamics of the physical particles. Consider the 
identity 

E «(- - «?> + E «(- - « = / |s«(« - , ) ' ilog(1 - eiZW) . (72) 

roots holes 7 

where the contour 7 is depicted in figure 6. Notice that the contour must not encircle poles of the 
function ~H(u), and therefore < e < 1. Then, we rewrite this contour integral as 

/°° flu f°° dv 

^-H'{u - v)Z(v) + / t^H'(u -v + ie) log(e- iZ ^ - 1) 
00 2?r J-00 27 ™ 



/CO J 
— H'(u - v - ie) log(l - e iZ ^) 
-00 27 ™ 



(73) 



We now use this result to substitute the sum over the Bethe roots in (72) by a sum over the holes, 
obtaining 



(l + —n')*Z(u) = -V h (u)+ nu-6 k )-^U(u-u k )-^2v(u-v k ) 

k=l k=l k=l 

_ -Lft'h*] * log(l - e iZ (»+ ie )) + — * log(e- iZ ( t '- ie ) - 1) . 
2ir 2tt 

1 17/ 



(74) 



Finally, after inverting the kernel 1 + j we obtain the final DdV equation 

^C") = ■7r'Peff(u) + ^ Meffiu-Ok) -^2V e ff(u-V k ) -^2U e ff(u-U k ) 

k=l k=l k=l (75) 

+ itfg/ * log(e- iZ ( u - ie ) - 1) - i/f^f 1 * log(l - e lZi - v+i ^) , 
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where 



K e ff{u) 
Veff(u) 

U eff (u) 



-i 



-i 



H 



eff 



-1 





:«)) 








:«)) 


F{V{u)) 






:«)) 








:«)) 






1 + T(±H> 


:«)) 



(76) 
(77) 
(78) 
(79) 
(80) 



where J 7 and J 7 1 denote the Fourier transform and inverse, respectively. The integrals in the 
Fourier transforms are again understood in the principal part sense. 

The convolutions appearing in this equation give the finite size corrections. For the moment, we 
will drop them, and study the infinite size limit of this equation. In the continuum limit, N — > 00 
and a — > 00, we define 

mL = 2Ne~™ (81) 
and therefore the source term in equation (76) can be rewritten as 

mL sinli7ra. (82) 

After taking the limit s — > 00 in the kernels (77-80) and exponentiating, we obtain the following 
result for very large mL 



#holes 

- imL sinh 7r6,- I I o2/n n 

= 11 ~ °k 

k^j k 



n 



0j -Vk- 



N R 

n 



9j — Uk — 2 



(83) 



The auxiliary Asymptotic Bethe equations can be obtained similarly. Specifying (68) for n = 2s + 1, 
we have 



N 



M 



V u (u k ) = 2nQf +1 + J2U(uj - uf) 



N R 



Nl 



(84) 



Vk 



k=l 



k=l 



The sum over the 2s-strings can be converted in a sum over the holes, using DdV (ignoring the 
finite size corrections) 



M Wholes 

Y J U{u-ul s ) = - U{u-d k ) + —U'*Z(;u 



(85) 



k=l 



k=l 
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Using the expression (75), we obtain after exponentiation 

fc 2 1 TT / 



ii; .. . i ii 



A similar equation could be derived for the strings of size 2s — 1. Therefore we have obtained the 
Asymptotic Bethe ansatz equations for the SU(2) PCF. 

We now turn to the second goal of this appendix, and consider the finite size corrections. We 
will restrict ourselves to the ground-state, in order to compare with the ground-state GKV equation 
(21). Hence, in the continuum limit and s — > oo, equation (75) can be rewritten as 

g 2 = -e imisinh ™exp (2ilm [K~ * log((#+) 2 - 1)]) . (87) 

where g 2 = e iZ ^ u \ and the shift in the contour of integration, e, was taken to be ^. This equation 
resembles the GKV equation (21), with the significant difference that denominator |g + | 2 — 1 is absent 
here. We assumed that the corrections to this hypothesis should be exponentially suppressed in 
N, and hence, negligible. But, as the number of Bethe roots is M = Ns, and also goes to infinity, 
the corrections are significant. Our assumption is then manifestly wrong. Therefore, this extra 
denominator can be interpreted as the string hypothesis deviation in the anti-ferromagnetic vacuum 
of the associated discretization (see section 5). 



B Details of the derivation of the lattice NLIE 

In this appendix, we give further details on the derivation of the hybrid equations (59). It is 
convenient to introduce the notation 

dlf = /d^lM e -*«. (88) 
J du 

With this notation, the contour integral is rewritten as 

dJT^e - * - dlT$,e* = . (89) 

Using (55) and (89) we get 



d/Q[±2*] = ®(^) ( ± dUl + ao ) e l T dl(l + a )e~^ + V dl^ ±2m± A , (90) 
2cosh l V ^ 1 

where O is the Heaviside step function. Next, we introduce a function clq ee ^ which is ANZC in 
the region |Im u\ <\. Using (54) we have 

2s-l 

dla = dlT 2s -i + dlQ [2s] e~ k - dlQ^ - ^ dl<f> [2n+2] . (91) 

n=0 

From (37) with j = 2s — 1, we express T2 S -i in terms of Y2 S _i 

1 1 2s ~ l 

tflTV-i = — — T dl(l + Y 2s _!) + — — T V dl(^-i +2 M- 2p+J - 2s] ) ■ (92) 
2 cosh % 2 cosh % ^— f 

2. 2. p=l 

We then plug (90) and (92) in (91). Finally, we invert the Fourier transform and integrate in u to 
get (59). 



2s-l 



22 



References 



[1] M. Luscher, "Volume Dependence Of The Energy Spectrum In Massive Quantum Field The- 
ories. 1. Stable Particle States," Commun. Math. Phys. 104, 177 (1986). 

[2] M. Luscher, "Volume Dependence of the Energy Spectrum in Massive Quantum Field Theories. 
2. Scattering States," Commun. Math. Phys. 105 (1986) 153. 

[3] A. B. Zamolodchikov, "Thermodynamic Bethe Ansatz In Relativistic Models. Scaling Three 
State Potts And Lee-yang Models," Nucl. Phys. B342, 695-720 (1990). 

[4] P. Dorey, R. Tateo, "Excited states by analytic continuation of TBA equations," Nucl. Phys. 
B482, 639-659 (1996). [hep-th/9607167]. 

[5] V. V. Bazhanov, S. L. Lukyanov, A. B. Zamolodchikov, "Integrable structure of conformal 
field theory, quantum KdV theory and thermodynamic Bethe ansatz," Commun. Math. Phys. 
177, 381-398 (1996). [hep-th/9412229]. 

[6] A. M. Polyakov and P. B. Wiegmann, "Theory of nonabelian Goldstone bosons in two dimen- 
sions," Phys. Lett. B 131, 121 (1983). 

[7] A. M. Polyakov and P. B. Wiegmann, "Goldstone Fields In Two-Dimensions With Multivalued 
Actions," Phys. Lett. B 141, 223 (1984). 

[8] L. D. Faddeev, N. Y. .Reshetikhin, "Integrability Of The Principal Chiral Field Model In 
(l+l)-dimension," Annals Phys. 167, 227 (1986). 

[9] A. Hegedus, "Finite size effects in the SS model: Two component nonlinear integral equations," 
Nucl. Phys. B679, 545-567 (2004). [hep-th/0310051]. 

[10] C. Destri, H. J. de Vega, "Light Cone Lattices And The Exact Solution Of Chiral Fermion 
And Sigma Models," J. Phys. A A22, 1329 (1989). 

[11] C. Destri, H. J. de Vega, "On The Connection Between The Principal Chiral Model And The 
Multiflavor Chiral Gross-neveu Model," Phys. Lett. B201, 245 (1988). 

[12] N. Gromov, V. Kazakov and P. Vieira, "Finite Volume Spectrum of 2D Field Theories from 
Hirota Dynamics," JHEP 0912, 060 (2009) [arXiv:0812.5091 [hep-th]]. 

[13] M. Beccaria, G. Macorini, "A numerical test of the Y-system in the small size limit of the 
SU(2)x SU(2) Principal Chiral Model," [arXiv: 1009.3428 [hep-th]]. 

[14] C. Destri, H. J. de Vega, "New thermodynamic Bethe ansatz equations without strings," Phys. 
Rev. Lett. 69, 2313-2317 (1992). 

[15] C. Destri, H. J. De Vega, "Unified approach to thermodynamic Bethe Ansatz and finite size 
corrections for lattice models and field theories," Nucl. Phys. B438, 413-454 (1995). [hep- 
th/9407117]. 

[16] A. B. Zamolodchikov and A. B. Zamolodchikov, "Factorized S-matrices in two dimensions as 
the exact solutions of certain relativistic quantum field models," Annals Phys. 120, 253 (1979). 

[17] N. Gromov, V. Kazakov, K. Sakai and P. Vieira, "Strings as multi-particle states of quantum 
sigma-models," Nucl. Phys. B 764, 15 (2007) [arXiv:hep-th/0603043]. 



23 



[18] Takahashi, M. 1971, Progress of Theoretical Physics, 46, 401 

[19] Takahashi, M. 1999, Thermodynamics of One-Dimensional Solvable Models, by Minora Taka- 
hashi, pp. 266. ISBN 0521551439. Cambridge, UK: Cambridge University Press, March 1999., 

[20] A. Zabrodin, "Discrete Hirota's equation in quantum integrable models," [hep-th/9610039]. 

[21] D. J. Gross and A. Neveu, "Dynamical Symmetry Breaking In Asymptotically Free Field 
Theories," Phys. Rev. D 10, 3235 (1974). 

[22] L. D. Faddeev, E. K. Sklyanin and L. A. Takhtajan, "The Quantum Inverse Problem Method. 
1," Theor. Math. Phys. 40, 688 (1980) [Teor. Mat. Fiz. 40, 194 (1979)]. 

[23] L. A. Takhtajan and L. D. Faddeev, "The Quantum method of the inverse problem and the 
Heisenberg XYZ model," Russ. Math. Surveys 34, 11 (1979) [Usp. Mat. Nauk 34, 13 (1979)]. 

[24] E. K. Sklyanin, "Quantum version of the method of inverse scattering problem," J. Sov. Math. 
19, 1546 (1982) [Zap. Nauchn. Semin. 95, 55 (1980)]. 

[25] N. Dorey and B. Vicedo, "A symplectic structure for string theory on integrable backgrounds," 
JHEP 0703, 045 (2007) [arXiv:hep-th/0606287]. 

[26] J. M. Maillet, "Hamiltonian Structures For Integrable Classical Theories From Graded Kac- 
Moody Algebras," Phys. Lett. B 167, 401 (1986). 

[27] J. M. Maillet, "New Integrable Canonical Structures In Two-Dimensional Models," Nucl. Phys. 
B 269, 54 (1986). 

[28] J. M. Maillet, "Kac-Moody Algebra And Extended Yang-Baxter Relations In The O(N) Non- 
linear Sigma Model," Phys. Lett. B 162, 137 (1985). 

[29] L. A. Takhtajan, "The picture of low-lying excitations in the isotropic Heisenberg chain of 
arbitrary spins," Phys. Lett. A87, 479-482 (1982). 

[30] H. M. Babujian, "Exact Solution Of The Isotropic Heisenberg Chain With Arbitrary Spins: 
Thermodynamics Of The Model," Nucl. Phys. B215, 317-336 (1983). 

[31] Zinn- Justin, P. 1998, eprint arXiv:solv-int/9810007, 10007 

[32] D. Volin, "Quantum integrability and functional equations," [arXiv:1003.4725 [hep-th]]. 

[33] Klumper, A., & Batchelor, M. T. 1990, Journal of Physics A Mathematical General, 23, L189 

[34] Klumper, A., Batchelor, M. T., & Pearce, P. A. 1991, Journal of Physics A Mathematical 
General, 24, 3111 

[35] P. P. Kulish, N. Y. .Reshetikhin, E. K. Sklyanin, "Yang-Baxter Equation and Representation 
Theory. 1.," Lett. Math. Phys. 5, 393-403 (1981). 

[36] P. P. Kulish, E. K. Sklyanin, "On the solution of the Yang-Baxter equation," J. Sov. Math. 
19, 1596-1620 (1982). 

[37] A. Klumper and P. Pearce, Conformal weights of RSOS lattice models and their fusion hier- 
archies, Physica A183 (1992) 304-350. 



21 



[38] A. Kuniba, T. Nakanishi and J. Suzuki, "Functional relations in solvable lattice models. 1: 
Functional relations and representation theory," Int. J. Mod. Phys. A 9, 5215 (1994) [arXiv:hep- 
th/9309137]. 

[39] A. Kuniba, T. Nakanishi and J. Suzuki, "Functional relations in solvable lattice model 2: 
Applications," Int. J. Mod. Phys. A 9, 5267 (1994) [arXiv:hep-th/9310060]. 

[40] J. Suzuki, "Spinons in magnetic chains of arbitrary spins at finite temperature," J. Phys. A 
A32, 2341-2359 (1999). 

[41] H. J. de Vega, F. Woynarovich, "Solution Of The Bethe Ansatz Equations With Complex 
Roots For Finite Size: The Spin S >= 1 Isotropic And Anisotropic Chains," J. Phys. A A23, 
1613 (1990). 

[42] A. Hegedus, "Nonlinear integral equations for finite volume excited state energies of the 0(3) 
and 0(4) nonlinear sigma-models," J. Phys. A A38, 5345-5358 (2005). [hep-th/0412125]. 

[43] N. Beisert and M. Staudacher, "Long-range PSU(2,2|4) Bethe ansaetze for gauge theory and 
strings," Nucl. Phys. B 727 (2005) 1 [arXiv:hep-th/0504190]. 

[44] N. Beisert, B. Eden and M. Staudacher, "Transcendentality and crossing," J. Stat. Mech. 0701 

(2007) P021 [arXiv:hep-th/0610251]. 

[45] N. Gromov, V. Kazakov and P. Vieira, "Exact Spectrum of Anomalous Dimensions of 
Planar N=4 Supersymmetric Yang-Mills Theory," Phys. Rev. Lett. 103 (2009) 131601 
[arXiv:0901.3753 [hep-th]]. 

[46] D. Bombardelli, D. Fioravanti and R. Tateo, "Thermodynamic Bethe Ansatz for planar 
AdS/CFT: a proposal," J. Phys. A 42 (2009) 375401 [arXiv:0902.3930 [hep-th]]. 

[47] N. Gromov, V. Kazakov, A. Kozak and P. Vieira, "Exact Spectrum of Anomalous Dimensions 
of Planar N = 4 Supersymmetric Yang-Mills Theory: TBA and excited states," Lett. Math. 
Phys. 91 (2010) 265 [arXiv:0902.4458 [hep-th]]. 

[48] G. Arutyunov and S. Frolov, "Thermodynamic Bethe Ansatz for the AdS§ x S 5 Mirror Model," 
JHEP 0905 (2009) 068 [arXiv:0903.0141 [hep-th]]. 

[49] A. V. Belitsky, "Fusion hierarchies for N = 4 superYang-Mills theory," Nucl. Phys. B 803 

(2008) 171 [arXiv:0803.2035 [hep-th]]. 

[50] A. Hegedus, "Discrete Hirota dynamics for AdS/CFT," Nucl. Phys. B 825 (2010) 341 
[arXiv: 0906. 2546 [hep-th]]. 

[51] N. Gromov, V. Kazakov, S. Leurent and Z. Tsuboi, "Wronskian Solution for AdS/CFT Y- 
system," arXiv: 1010.2720 [hep-th]. 

[52] V. Kazakov and S. Leurent, "Finite Size Spectrum of SU(N) Principal Chiral Field from Dis- 
crete Hirota Dynamics," arXiv: 1007. 1770 [hep-th]. 



25 



